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Abstract: Several models for physics beyond the Standard Model predict new particles 
with a decay signature including hard photons and missing energy. Two well-motivated 
examples are supersymmetry with gauge-mediated breaking (GMSB) and the standard 
model with two universal extra dimensions. Both models lead to decay chains with sim- 
ilar collider signatures, including hard photon emission. The main discriminating feature 
are the spins of the new particles. In this paper we discuss how information about the 
spins of the particles can be extracted from lepton-photon or quark-photon invariant mass 
distributions at the Large Hadron Collider. The characteristic shapes of the distributions 
are derived analytically and then studied in a realistic Monte-Carlo simulation. We find 
that for a typical GMSB mass spectrum with particle masses below 1 TeV, already 10 fb~^ 
integrated luminosity at 14 TeV center-of-mass energy are sufficient to discriminate the 
two models with high significance. 
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1. Introduction 

Several models for physics beyond the Standard Model (SM) introduce partner particles 
for all SM particles, the lightest of which is stable and could be the constituent of dark 
matter. The best known examples are supersymmetry (SUSY) with conserved R-parity 
and universal extra dimensions (UED) with conserved Kaluza-Klein (KK) parity ||^] . Since 
the observable signatures for these models at the Large Hadron Collider (LHC) look quite 
similar it will be important to test the fundamental quantum numbers of the new 
particles in order to scrutinize the nature of the underlying physics. 

For example, both SUSY and UED require that the couplings of the new particles 
are identical to the corresponding couplings of their SM partners; a prediction which can 
be tested at the LHC by measuring cross section ratios [||. However, a crucial distinction 
between the two models is given by the spins of the new particles. While the SUSY partners 
differ from their SM counterparts by one half-unit of spin, the KK excitations in UED 
have the same spin as their SM partners. Recently, extensive work has been performed to 
determine the spins of SUSY or UED particles by exploiting angular correlations in the 
decay of those particles at the LHC^. Many papers have focused on decay chains involving 
lepton pairs (sj, @, IC , 11,12]- A typical example of such a decay chain in SUSY is 



q ^ qX2 ^ ql^^ Xi- Other studies have examined channels involving heavy gauge bosons 
||9|, |l3|, 14, 15|, sleptons |14, 16] and top quarks. However, in all of the existing studies it 



^ A few studies have explored more model-dependent discrimination methods based on total cross sections 
and higher KK modes M. 
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was assumed that the hghtest new particle is either a neutrahno or weak vector boson in 
SUSY or UED, respectively. 

In this article the spin determination from angular correlations is extended to decay 
chains that involve hard photons in the final state. Such decay channels occur naturally 
in gauge- mediated supersymmetry breaking (GMSB), where the lightest SUSY particle is 
the gravitino as well as in the extension of the SM by two universal extra dimen- 
sions (UED6), where the lightest KK mode is typically a scalar component of a higher- 



dimensional vector boson, called "scalar adjoint" p9| , [20| , 21 1. It has been shown earlier 
that a high-energy e+e" collider could distinguish between GMSB and UED6 by studying 
angular correlations in pair production and decay of the selectron (KK-electrons) and neu- 
tralinos (KK-gauge bosons) [^. The purpose of the present paper is to study how such 
a distinction can be achieved at the LHC by analyzing decay chains involving leptons and 
photons. In particular, we are investigating how the spin of the lightest new particle, which 
escapes from the detector in form of missing momentum, can be inferred from invariant 
mass distributions of the leptons and photons. 

In section ^ we describe analytical calculations of the relevant invariant mass distri- 
butions and compare the angular correlations predicted by GMSB and UED6. In order to 
evaluate the prospects for experimental measurements of these distributions, we present in 
section |3| results of a realistic Monte-Carlo simulation, incorporating the spin correlation 
effects of the two models. We present our conclusions in section ^ The detailed analytical 
results are collected in the appendix. 

2. Spin correlations in GMSB and UED6 



In our notations and conventions we follow Ref. [^3[ for supersymmetry and Ref. 2C, 21 1 
for the Standard Model in six dimensions (UED6). 

GMSB is a promising candidate for a mechanism that generates TeV-scale masses for 
the SUSY partners, most notably since it explains the absence of large flavor-changing 
neutral currents. In GMSB, the lightest SUSY particle is typically the gravitino G, with 
a mass itIq < 1 MeV. If the next-to-lightest SUSY particle is a neutralino, it will mostly 
decay into the gravitino and a photon, Xi ~^ jG. Depending on the neutralino mass, there 
could be a smaller branching fraction into a Z boson, which we will not investigate further. 

In UED6, the lightest particle at KK level (1,0) is typically the scalar adjoint of the 
hypercharge boson, B^^'^^ [20, 24]. The vector mode of the KK-hypercharge boson, Bj^'^\ 



can decay into the scalar adjoint via a loop-induced process, B^^'^^ B^j^'^\ This decay 



mode has a sizable branching fraction of about 34% Therefore the two models lead 
to very similar decay signatures. In particular, the decay chain 

^o^/±[T^/+rx;^/+r7G, (2.1) 

which is typical in GMSB, is imitated by the equivalent process in UED6, 
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Figure 1: SUSY (red/top) and UED6 (blue/bottom) decay chains with observable Gnal state 

see Fig. |l|. Both processes lead to a final state signature of a same-flavor, opposite-sign 
lepton pair, one photon, and missing transverse momentum ^-p. 

In the following we will study spin correlation effects in these decay chains in detail. 
We will also briefly analyze the following shorter decay modes of squarks: 

QR^qxI^q-fG (2.3) 

for GMSB, and 

Qa,o)_ ^^(1,0)^^^^(1.0) (2.4) 

for UED6. Although these decay channels have considerably larger SM backgrounds at the 
LHC, they also have larger branching ratios compared to the decay channels with leptons. 
They could be very useful to obtain information on the spin of the quark partners. 

While both GMSB and UED6 have similar mass hierarchies, which allow the decay 
channels in eqs. ( |2.lD -(^^), the typical mass spectra are quantitatively rather different. 

In minimal gauge mediation the gaugino mass parameter relation Mi/gf = = 
M^/g'^ implies that the weak gauginos are rather light, while the gluino is much heavier. 
Furthermore, the squarks are also heavy, while the sleptons have masses close the the 
gauginos. For our numerical analysis in the next section we will use the reference scenario 
Gla from Ref. ||2^. The masses of the particles appearing in our decay chain and their 
branching ratios are summarized in Tab. |l[ 

In universal extra dimensions, on the other hand, the masses of all particles of one 
KK level (j, k) have the same value = (j^ + k'^)/B? , where R is the size of the extra 
dimensions. This degeneracy is lifted only by radiative corrections |20, 24], which shift the 



masses by up to 20%. As a result, a typical spectrum for the particles of KK level 1 in 
UED6 is much less hierarchical than in GMSB. In our numerical analysis, we will use the 
reference scenario Ul, defined by = 500 GeV, with the masses and branching ratios 
|2l|] given in Tab. |l|. 

These qualitative features of the spectra could be used to distinguish GMSB and UED6 
experimentally. However, there are several caveats to consider: The mass spectra of GMSB 



can vary substantially in non-minimal models, see e. g. Ref. |26|. Extra dimensional models 
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Table 1: Masses and branching fractions for the GMSB scenario Gla (left) and the UED6 scenario 
Ul (right). 



are known to become strongly coupled at large energies and require some new physics to 
be present at the scale. The effects of this unknown high-scale physics could generate mass 



contributions to the KK particles in UED |27|. However, the spins of the the new particles 
can serve as very robust discriminators between different models. 

A non-zero spin of a particle can lead to angular correlations between its decay prod- 
ucts. At the LHC, angular correlations are manifested in the invariant mass distributions 
of the visible decay products of a decay chain. The long decay chains (2.1) and ( |2.2D are 
of the general form 



itljB 



(2.5) 



with uiA < ms < rric < mo. Here we call the lepton that is emitted in the first decay step 
the "near-lepton" while the lepton from the second decay step is named the "far-lepton" 
Prom this one can construct the invariant masses 
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(C) 

which are related to 0„ / , the angle between the near-lepton and the far-lepton in the rest 

{ H) 

frame of C, 9n'y , the angle between the near-lepton and the photon in the B rest frame, 
( R") 

and Gf^ , the angle between the far-lepton and the photon in the B rest frame, respectively. 
The maximum values for the invariant masses are given in eq. (A. 2) in the appendix. 



For a given decay matrix element, the distribution with respect to some invariant mass 
is then obtained by integrating over all remaining kinematical variables in a given reference 
frame, as described in detail in Refs. [0, 28|. 



In practice, the near and far leptons cannot be distinguished in a straightforward way. 
Instead the observable lepton-photon invariant mass distribution is the sum of dT / dm^^ + 
dr/dm^^. Additional information can be obtained from the distributions with respect to 
the minimum and maximum of the lepton-photon invariant masses, 

ml^ = max{m^^, mj^}, mf^ = min{m^^, mj^}. (2-10) 

Since the total magnitude of the decay width does not carry any information about the 
spins of the particles involved, we will normalize the invariant mass distributions to unity, 

Tq dm dm ' ^ ^ 

where Tq is the integrated decay width of the given decay channel, and dP is defined as a 
differential probability density. 

For the short decay chains (^]^) and ( |2.4D of the general form 

C ^ qB ^ q'jA, (2.12) 

with mA < ms < nic, the only observable invariant mass distribution that can be con- 
structed is 

2 / . n2 / rna,x\2 1 



Here {m^^'^)'^ = {m^ — m^)(m^ — m\)/m'^. 
2.1 Spin correlations in GMSB 

In GMSB the final state l~^l~'y + ^-f is fed by the decay chain eq. ( |2.lD , with and Xi &s 
intermediate particles. Since the is a scalar it does not transmit any angular correlations. 
The fermionic x? can lead to non-trivial spin correlations in the decay chain, but only if the 
couplings at both the production and decay vertex are chiral, i. e. left- and right-handed 
components have different coupling strength Q. While this condition is fulfilled for the 
l^-P-Xi vertex, the XrJ-G vertex has the form ||2^ 




oc 



^',7M]7a. (2.14) 
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(VFVS) 



(SFVS) 




Figure 2: Different UED6 decay chains witfi various spin configurations leading to tie same 
+ ^j- signature. 



Here the left- and right-handed component of the neutrahno couple with equal strength 
so that no observable angular correlations are generated. Thus the only features in the 



invariant mass distributions of the GMSB decay chain {2A) are generated by the phase 
space. Analytical results for the lepton-photon and lepton-lepton distributions are listed 
in the appendix. 

Similarly, the short decay chain eq. (|2.3| ) does not lead to visible angular correlation 
effects. The quark-photon invariant mass distribution is given by 

dP _ m| 
dm2^ (ml - m%){ml - m\) ' 



By comparison with the formulas in the appendix one can see that eq. ( p. 15 ) is identical 
to the rrify distribution of the long chain. This can be easily understood by the fact that 
the chirality of the slepton and squark couplings is identical. 



2.2 Spin correlations in UED6 

The typical mass hierarchy generated by radiative corrections in UED6, m (i,o) > m (i,o) > 



m (1,0) > m (1,0) , enables the decay chain eq. (|2.2|) . However, in general sizeable corrections 
to the KK-particle masses could be generated by the unknown physics that complete the 
theory at high energies [27|. Thus for completeness we will study all possible decay chains 
that, for suitable mass hierarchies, could lead to the final state + ^j, as illustrated 

in Fig. |, 
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^ /± ^ /+/- ^i-o) ^ /+r 7i?(;'°^ (VFVS), (2.16) 

4!'°) ^ /± - 4^'°) ^ /+r 7i3;j'°^ (SFVS), (2.17) 

^ /± L+p(i'0) ^ /+r ^ l+r 7 ^I'O) (VFSV), (2.18) 

^I'O) ^ /± ^ /+/- ^ 7i3(^'°) (SFSV). (2.19) 

Here we have introduced short-hand notations for the four decay chains based on the KK 
particles at each decay stage being a scalar (S), fermion (F) or vector (V). In all cases 
we keep the couplings of the KK particles as they are predicted by the UED6 model, and 
mixing between gauge eigenstates is neglected. Analytical results for the invariant mass 
distributions for all four combinations are listed in the appendix. 

In case of the short decay chain, eq. ( |2.4| ), there are two possible decay chains with the 
same final state, depending on the mass hierarchy. 



q(1,0) ^ ^41,0) ^ ^^^(1,0) (psy)_ (2.21) 



The quark-photon invariant mass distributions for the two cases read 



FVS : 

dP 3m^ (2m^^m| - 2m^^(m^ - m|)(m^ - m^) + {m\ - m'^^)'^{m^ - m^)) 
dm^^ (m^ — m^)3(m^ — m^)2(2m^ + m^) ' 

(2.22) 



FSV : 

dP m 



B 



dm^^ (m^ — m\)(jn\ — m\) 



(2.23) 



As before, one can see that these are identical to the rnj^ distributions of the long chain for 
the VFVS/SFVS and VFSV/SFSV combinations, respectively. In the FSV case, since the 
intermediate B particle is a scalar, this decay process does not involve any spin correlations 
and is identical to the pure phase space distribution, and thus to the GMSB distribution. 

2.3 Discussion of analytical results 

Figs. ^ and ^ show the distributions for the four independent observable invariant mass 
combinations of the final state: the di-lepton invariant mass m„j, the "low" and 

"high" lepton- invariant masses m;^ and m/i^, respectively, and the lepton-lepton-photon 
invariant mass 7ti„j^. Each plot contains five curves corresponding to the five models 
(or spin assignments) GMSB, VFVS, SFVS, VFSV, and SFSV. In case of Fig. |, for all 
five models the masses have been chosen from the Gla scenario, with ttia = rn^ = 0, 
TUB = iT^^o, mc = mg^, and mu = m^o. On the other hand, Fig. ^ shows the situation for 

the Ul spectrum with niA = m^nm, tub = m^nm, mc = m^}'^\ and nif) = m^ri.o). 

i>jj ^+ 

As evident from the plots, different distributions could discriminate between different 
spin assignments. The di-lepton distribution dP/dm^j is markedly different for the VFSV 
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Figure 3: Observable invariant mass distributions for the decay chain in eq. (2^) for different 
models and masses from the Gla scenario. 



and SFSV models, where the B particle is a scalar, compared to the other models. This 
can be understood from the fact that the chiral structure of the KK-fermion couplings to 



KK-scalars or KK- vector bosons is different (see for example Ref. |10| for more details). 
On the other hand, the peak in the dP/dm/j^ and dP/drUn/'y spectra is relatively enhanced 
for the VFSV and SFVS models. As a result, the discriminative power between different 
spin assignments is maximized by including all four distributions in the analysis. 

Fig. ^ shows the jet-photon invariant mass distribution for the short decay chain. For 
the distribution on the left hand side the masses have been chosen from the Gla scenario, 



rriA 



m 



G 



0, niB 



mc 



and on the right hand side the distribution 



is shown for the Ul scenario, where niA = "^r(i,o), ms = m (i,o), mc = m (i.o). As 
mentioned before, the GMSB and FSV chains both do not generate any spin correlations 
and thus cannot be distinguished from each other. However, the plots show a discrimination 
potential between FVS and the other two cases, especially for the Gla mass scenario. The 
qualitative features of the jet-photon distribution are similar for the Gla and Ul scenarios, 
although the spin correlation effects are less pronounced for the more degenerate Ul mass 
spectrum than for the more hierarchical Gla spectrum. 



3. Monte-Carlo simulation and numerical analysis 

In the previous section theoretical formulas for the invariant mass spectra of decay chains 
in GMSB and UED6 were discussed. However, in a realistic experimental setup one must 
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Figure 4: Observable invariant mass distributions for the decay chain in eq. (2.5) for different 
models and masses from the Ul scenario. 



take into account several effects that influence the measured distributions. 

• Detector efi^ects, such as hmited resolution, and reconstruction efi^ects from particle 
identification can smear out the visible invariant mass distributions. Further, detector 
acceptance can distort the shape of the distribution. 

• For the separation of the new physics signal from the SM backgrounds suitable selec- 
tion cuts have to be implemented, which can distort the shape of the distributions. 
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Figure 5: Observable invariant mass distributions for the short decay chain in eq. ( 2.12| ) for masses 
from the Gla scenario (left) and Ul scenario (right). 
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• Since both SUSY and KK particles are produced in pairs, often there wih be two 
hard photons in one event, leading to a two- fold ambiguity in the reconstruction of 
the correct decay chain. 

Therefore we have performed a realistic experimental simulation of the new physics signal 
from the long decay chain for GMSB and UED6, with the goal of reconstructing the 
spin-sensitive invariant mass distributions and discriminating between models from the 
simulated data. 

The decay chains in eqs. (I2II ), (|2T^ )- (|2lg| ) have been computed with CompHEP 4.4 



|30f| for 14 TeV proton-proton collisions, using the GMSB model file from Ref. |31] and our 



own implementation of the UED6 model. Parton-level events generated with CompHEP 



were then passed on to Pythia 6.4.12 |32|. Finally, the ATLAS detector] 33] is simulated 
using the parametrized fast detector simulation ATLFAST[|^, which includes detector 
acceptance, resolution and some basic particle identification. 

For the Gla and Ul mass spectra, the aforementioned decay chains are initiated mostly 
by squark (KK-quark) and gluino (KK-gluon) production processes, respectively. Since 
our analysis does not depend crucially on the details of the hadronic decay products of the 
squarks (KK-quarks) and gluinos (KK-gluons) we have only generated events for squark 
(KK-quark) pair production as the primary hard process, and then normalized the total 
event count according to the total production cross section including gluino (KK-gluon) 
processes. 

As a further simplification, we have only generated matrix elements with the decay 
chain of one of the two squarks (KK-quarks) within CompHEP, which correctly implements 
all spin correlations. The decay of the second squark (KK-quark) was simulated in Pythia, 
without spin correlations. Nevertheless, this procedure provides a good approximation to 
the complete matrix elements since the branching ratio for the decay chain leading to the 
final state I'^l"^ is relatively small and thus only a very small fraction of the events contains 
two decay cascades of this type. Moreover, Pythia was used for simulating initial and final 
state radiation and hadronization. 

The event selection was performed according to Ref. [25|. First the effective mass 

MeS = ^T+PT,l+PT,2+PT,3+PT,4 (3.1) 

is defined, where ^-f- denotes the missing transverse energy and pT,i the transverse momenta 
of the 4 hardest jets. The selected events need to fulfill the following conditions: 

1. 4 jets with transverse momenta > 25 GeV, 

2. Meff > 400 GeV, 

3. >0.1Meff, 

4. 2 photons with transverse momenta > 20 GeV, 

5. 2 electrons or muons with transverse momenta > 20 GeV. 
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After application of these cuts the SM background is reduced to a neghgible level |]25| , 
while about 20% of the signal is retained. 

Since the selected events include two photons, one of them has to be selected to compute 
the invariant mass distributions. Good results are obtained, when choosing the one photon, 
which gives the smaller iT^nf-y invariant mass. 

It would be interesting to also analyze the short decay chain in a complete simulation, 
since the expected rates are large (about 5 pb both for Gla and Ul). However, the signature 
of this final state, two hard jets, two photons, and missing energy, is very sensitive to issues 
related to jets faking photons. Therefore, it would require a more careful analysis of QCD 
backgrounds, which will be left for a future publication. 

3.1 Cross sections 



For a realistic analysis, the model-dependent cross-sections and event numbers have to 
be calculated. At the Gla point, the total squark and gluino production cross section is 
aq/g = 7.6 pb ^5|. Since the gluino is lighter than the squarks, the squarks decay through 
cascades involving a gluino. With the branching ratios from Tab. |l] a good approximation 
for the cross section for the decay chain is 



tTGia = <Tq/g X 2 X 2 X BR ^ qqxl] x BR [xl ^ e+e'x?] x BR 
~ 1.2 pb. 



X? 



(3.2) 



The factors of 2 result from the two generations of leptons in the di-leptonic decay chain 
and the fact that squarks or gluinos are produced in pairs and both of them can decay 
through this channel. 

In the Ul model, the KK-quarks and KK-gluons are produced with the following cross 



sections 21 



O" (1,0) „(i,o) 



fT^(l,0) (1,0) ~ 



7pb, 
10 pb, 
26 pb. 



18 pb, 
24 pb. 



(3.3) 



In this scenario, KK-quarks are lighter than the KK-gluon, such that the latter will mostly 
decay into KK-quarks. Therewith, the total KK-quark production cross section is 



^■^(i.o) — 2 0-^(1, o)q(i.o) + '^q(i.o)q{i>o) + 2 BR 



a (1,0) (1,0) 



(3.4) 



+ 1 + BR 



O"^(1,0)q(1,0) + BR 



(7(1,0) 



fT^(l,0)Q(l,0). 



As above, factors of 2 account for the two sides of the pair production process. Then the 
total cross section for the decay chain adds up to 



cj (1,0) X 2 X 2 X BR 



Q 



(1,0) 



X BR 



r(i,o) 



0.12 pb, 



X BR 



^(1,0) 



(3.5) 
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Vr V b 


Vr o V 


or Vo 


or o V 
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0.000 (m,,^) 


0.000 (m„/) 


0.006 (mh^) 


0.000 (mnf) 


VFVS 


0.056 {rrih^) 




0.000 (ninf) 


0.000 (m/,^) 


0.000 


VFSV 


0.577 (m„/) 


0.155 (nih'y) 




0.000 (mn/) 


0.000 (m„/) 


SFVS 
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0.065 (m,^) 


0.084 (m/i^) 




0.000 (m„/) 


SFSV 


0.000 (m„/-) 


0.000 (m„/) 


0.000 (m„/) 


0.000 {rrinf) 





Table 2: Minimal -probabilities for 10 fb^^. Tiic distributions that provide the strongest con- 
straints are noted in parentheses. The values in the blue cells (upper right) are for the Gla mass 
spectrum and in the red cells (lower left) are the values for the Ul mass spectrum. 



Using the signal efficiency of the cuts in Ref. |25] and assuming an integrated luminosity of 



10 fb~^ we obtain NqiI = 2500 selected events for the Gla point and N^^^ = 250 selected 
events for the Ul model. 

It should be noted that the event rates depend strongly on the underlying model and 
its parameters, as well as the choice of selection cuts. Furthermore, the signal efficiency 
after cuts depends also on the decay chain of the second squark or KK-quark, whose 
branching ratios vary between different models. For the purpose of this study, we do 
not vary the choice of cuts, cross sections and signal efficiency between different models 
when we compare them for one given mass spectrum (i. e. for all spin assignments we 
assume 2500 selected events for the Gla mass spectrum and 250 selected events for the Ul 
mass spectrum). Rather, our numerical analysis of the two scenarios should only serve as 
concrete examples for a spin determination of a new physics signal, in particular since our 
method does not rely on information about total rates. 

3.2 analysis 

In order to discriminate the histograms for the different spin configurations we used the x^- 
test implemented in ROOT |36|. It returns the x^-probability, i. e. the probability that two 



histograms with identical underlying distribution functions have a bigger value than the 
two compared ones. Since these values depend on the number of bins, the discrimination 
was performed with 5 bins, which showed the best discriminative power at 10 fb-^ In 
Tab. |2| the minimal x"^ probabilities for each pair of spin configurations are listed. 

The results of the x"^ probabilities reflect the general features that can be seen in 
the histograms of Fig. |6| and 0. As expected, the discrimination between different spin 
combinations is far more effective for the Gla scenario than for the Ul scenario, as a result 
of the more degenerate mass spectrum and lower cross section in the latter case. For the 
Gla scenario, even with 10 fb~^ almost all models can be distinguished with a confidence 
level better than 99.9%. In the case of the Ul scenario only the SFSV spin assignment can 
be distinguished at this confidence level with 10 fb~^ luminosity, owing to the distinctly 
different shape of the distribution of the di-lepton invariant mass In other cases 

however, in particular for the GMSB spin assignment and the extra-dimensional chain 
VFSV, it is not possible to make a distinction even at the 95% confidence level. Here a 
much larger integrated luminosity would be required for a significant discrimination. 
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Figure 6: Reconstructed invariant mass distribution from ATLFAST detector simulation, corre- 
sponding to Gla masses and cross sections for 10 fb~^ integrated luminosity. The histograms have 
been divided into 5 bins. 



As an example Fig. |8| shows the reconstructed invariant mass distributions correspond- 
ing to Gla masses for an integrated luminosity of 30 fb~^. The distributions are divided in 
10 bins, which gives the best discriminative power for this luminosity. More shape informa- 
tion is available due to more bins and events allowing for a better separation between the 
different models. For the Ul mass scenario the GMSB spin assignment can be distinguished 
from the extra-dimensional chain VFSV at almost 90% confidence level. 

Nevertheless, in summary the analysis of invariant mass distributions proves to be a 
powerful tool to identify the spins of particles in decay chains with hard photons. In this 
analysis the form the of the heavy particle interactions was fixed to the GMSB and UED6 
predicitions. In principle, the shapes of the invariant mass distributions are also sensitive 
to the couplings of those particles [^] , but a more general analysis, where these parameters 
are taken as free variables, will be left for future work. 
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Figure 7: Same as Fig. but for Ul masses and cross sections. 



4. Conclusion 



In this article we have analyzed the prospects for determining the spins of new particles in 
decay chains with photons and missing energy at the LHC. As concrete model incarnations 
of such signatures we considered supersymmetry with gauge-mediated breaking (GMSB) 
and the standard model with two universal extra dimensions (UED6). Each of these models 
predicts new partners of the SM particles, the lightest of which is stable on grounds of a 
conserved parity. The super symmetric or KK partner of the photon can decay into this 
stable particle and a hard photon. At the LHC the partners of the colored particles 
(squarks/gluinos or KK-quarks/KK-gluons) are produced with large cross sections and 
subsequently could decay in several steps until the final photon emission step. As became 
well known in recent years, the distributions of the invariant masses of two or more of the 
decay products is sensitive to the spins of the decaying particles. The measurement of the 
spins of the new particles is of central interest since this is the crucial difference between 
supersymmetry and extra dimensions. 

Following this approach, we first derived analytical expressions for the invariant mass 
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Figure 8: Reconstructed invariant mass distribution from ATLFAST detector simulation, corre- 
sponding to Gla masses and cross sections for 30 fb~^ integrated luminosity. The histograms have 
been divided into 10 bins. 



distributions for characteristic decay chains of those two models. In greater detail we an- 
alyzed a class of decay chains that lead to two leptons, one photon, and missing energy. 
This signature can stem from five different decay processes (one in GMSB, four in UED6), 
which differ by the spins of the intermediate particles. It was found that different invariant 
mass distributions show different distinctive features between those five cases, so that the 
discriminative power is maximized by combining information from all distributions. In a 
second step, we performed a realistic phenomenological Monte-Carlo analysis including a 
fast detector simulation for those five processes. Due to detector effects and cuts for back- 
ground rejection, the reconstructed invariant mass distributions are distorted compared to 
the analytical parton- level results. Nevertheless, the essential characteristic features of the 
five different spin assignments are preserved, so that the spins of the intermediate particles 
can be determined. 

For a typical GMSB mass spectrum with superpartner masses below 1 TeV, we found 
that with 10 fb~^ integrated luminosity almost all of the five different spin assignments, 
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and in particular the GMSB and UED6 models, can be distinguished with 99.9% confidence 
level. A typical UED6 spectrum, on the other hand, is more degenerate, leading to smaller 
branching ratios to leptons and to suppressed spin correlation effects in the invariant mass 
distributions. As a result, for such a mass spectrum, GMSB and UED6 can only be distin- 
guished with a confidence level of less than 50% with 10 fb~^. A much higher luminosity 
would be needed for better discrimination in this scenario. 

Additional information could be obtained by looking at other decay chains. We briefly 
investigated decay chains leading to a hard jet, a photon, and missing energy, but no 
leptons. Such a decay chain has a large branching fraction for UED6 scenarios and thus 
might be useful for discriminating models in this case. A more conclusive answer would 
require a detailed analysis of backgrounds to this process, which is beyond the scope of 
this work. 

In summary, it appears feasible to distinguish GMSB and six-dimensional UED models 
with LHC data alone if the mass spectrum is not very degenerate. 
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A. Analytical results for invariant mass distributions 

In the appendix we list analytical results for the invariant mass distributions of the various 
GMSB and UED6 decay chains discussed in section H. All those decay chains are of the 
form 

D^ltC^ltljB^ltlJjA, (A.l) 

with rriA < ms < mc < mo. Owing to kinematical constraints, the invariant mass 
distributions are divided into sections, which are bounded by the kinematical edges 
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